UHE NEUTRINO ASTRONOMY AND NEUTRINO OSCILLATIONS 



V. Berezinsky 
INFN, Laboratori Nazionali del Gran Sasso, 
67010 Assergi (AQ) Italy 



ABSTRACT 

UHE neutrinos with E > 10 17 eV can be produced by ultra-high energy cosmic 
rays (UHECR) interacting with CMB photons (cosmogenic neutrinos) and by 
top-down sources, such as topological defects (TD), superheavy dark matter 
(SHDM) and mirror matter. Cosmogenic neutrinos are reliably predicted and 
their fluxes can be numerically evaluated using the observed flux of UHECR. 
The lower limit for the flux is obtained for the case of pure proton composition 
of the observed UHECR. The rigorous upper limit for cosmogenic neutrino flux 
also exists. The maximum neutrino energy is determined by maximum energy of 
acceleration, which at least for the shock acceleration is expected not to exceed 
10 21 — 10 22 eV. The top-down sources provide neutrino energies a few orders of 
magnitude higher, and this can be considered as a signature of these models. 
Oscillations play important role in UHE neutrino astronomy. At production of 
cosmogenic neutrinos r-neutrinos are absent and D e neutrinos are suppressed. 
These species, important for detection, appear in the observed fluxes due to 
oscillation. Mirror neutrinos cannot be observed directly, but due to oscillations 
to ordinary neutrinos they can provide the largest neutrino flux at the highest 
energies. 

1. Introduction 

The boundary E v ~ 10 17 eV between HE and UHE neutrino astronomy is con- 
nected with observational technique. The neutrino observations with underwater /ice 
detectors are valid mostly for HE neutrino range. This is successful and well de- 
veloped experimental technique. Two detectors, Baikal and AMANDA, has reached 
now the saturated statistics, two other arrays, IceCube and ANTARES, started the 
observations, and two projects, NESTOR and NEMO, are at the stage of testing. 

The prospects for UHE neutrino astronomy appeared in 1960s soon after predic- 
tion of the GZK D cutoff. It has been realized^ 1 that proton interaction with CMB 
photons at large redshifts in case of cosmological evolution of the sources can produce 
UHE neutrino fluxes much higher than the observed UHECR flux. 

In 80s it was understood that topological defects can produce unstable superheavy 
particles with masses up to the GUT scale ® and neutrinos with tremendous energies 
can emerge due to this process®. 

It has been proposed that UHE neutrinos can be detected observing the horizontal 
Extensive Air Showers (EAS) The exciting prospects for detection of UHE neutri- 
nos have appeared with the ideas of space detection, e.g. in the projects EUSO® and 
owiR At present there is well developed JEM-EUSO project ^ with the prospects 
to start the observations in 2012 - 2013. 



The basic idea of detection by EUSO is similar to the fluorescent technique for 
observations of extensive air showers (EAS) from the surface of the Earth. The UHE 
neutrino entering the Earth atmosphere produces an EAS. The known fraction of 
its energy, which reaches 90% , is radiated in form of isotropic fluorescent light. An 
optical telescope from the space detects it. Since the observatory is located at very 
large height (~ 400 km) in comparison with thickness of the atmosphere, the fraction 
of detected flux is known, and thus this is the calorimetric measurement (absorption 
of up-going photons is small). A telescope with diameter 2.5 m controls the area 
~ 10 5 km 2 and has a threshold for EAS detection E t h ~lx 10 19 eV^*. 

The very efficient method of UHE neutrino detection is given by observations of 
radio emission by neutrino-induced showers in ice and lunar regolith. This method 
has been originally suggested by G. Askaryan in 60s®. Propagating in the matter the 
shower acquires excessive negative electric charge due to involvement of the matter 
electrons in knock-on process. The coherent Cerenkov radiation of these electrons 
produces the radio pulse. Recently this method has been confirmed by the laboratory 
measurements There were the experiments to search for such radiation from 
neutrino-induced showers in the Greenland and Antarctic ice and in the lunar regolith. 
In all cases the radio-emission can be observed only for neutrinos of extremely high 
energies. The upper limits on the flux of these neutrinos have been obtained: in 
GLUE experiment by radiation from the moon, in FORTE experiment ^ by 
radiation from the Greenland ice and in ANITA and RICE 03 experiments from 
the Antarctic ice. 

The characteristic feature of these detection methods is the high energy threshold 
E > 1 X 10 19 — 1 x 10 20 eV. How neutrinos of these energies can be produced? 

The most conservative mechanism of UHE neutrino production is p7 mechanism 
of collisions of accelerated protons/nuclei with low-energy CMB photons. To provide 
neutrinos with energies higher that 1 x 10 20 eV the accelerated protons must have 
energies higher than 2x 10 21 eV. For non-relativistic shock acceleration this energy can 
reach optimistically 1 x 10 21 eV. For relativistic shock this energy can be somewhat 
higher. To have neutrinos with higher energies one has to put his hopes on less 
developed ideas of acceleration such as acceleration in strong e-m waves, exotic plasma 
mechanisms of acceleration and unipolar induction. 

The top-down scenarios can naturally provide neutrinos with energies higher and 
much higher than 1 x 10 20 eV. The idea common for many mechanisms is given by 
existence of superheavy particles with very large masses up to GUT scale. In Grand 
Unified Theories (GUT) these particles (gauge bosons and higgses) are short-lived. 
In the cosmic space they are produced by Topological Defects (TD). The decay of 
these particles results in the parton cascade, which is terminated by production of 
pions and other hadrons. Neutrinos are produced in their decays. 

The superheavy particles are naturally produced at post-inflationary stage of the 
universe. The most reliable mechanism of production is gravitational one. The masses 



of the produced particles can reach 10 13 — 10 14 GeV. Protecting by some symmetry 
(e.g. gauge symmetry or discrete gauge symmetry like R-parity in supersymmetry), 
these particles can survive until present time and produce neutrinos in the decays or 
annihilation. 



2. Effects of UHE neutrino oscillations 



Characteristic distances to the sources of UHE neutrinos are much larger than 
maximum oscillation lengths. Taken as the maximum neutrino energy -E max ~ 10 20 eV 
one obtains the maximum oscillation length of order of 100 pc 
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while the distances r to UHE neutrino sources even in our galaxy is of order of a few 
kpc. The relation r ^> £ osc is often interpreted as a sufficient condition for flavour 
equipartition at observation u e : : v T = 1 : 1 : 1 . However, it is well known that 
this is not true for arbitrary flavour composition at generation. In the recent work 
the precise calculations for connection of the generation and observed flavour 
compositions are performed. 

The basic idea of calculations in can be explained in the simplified form as 
follows. 

Let us consider how observed flavour ratio is connected with flavour composition 
at generation in the limiting case r ^> £ osc . The flavour neutrino eigenstate u a is 
given by mixing of mass eigenstates U{ as v a = U a kVk, where a = e, /x, r, are flavour 
indicies, k are mass eigenstate indicies k — 1,2, 3 and U a k is the mixing matrix. In 
the most general case U a k is expressed through solar #i 2 and atmospheric 623 neutrino 
mixing angles, and include also small term (see equation below). However with 
a good accuracy one may use U e s = and 823 = vr/4 and thus obtain 
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In the limiting case r ^> £ osc , < sinr// osc >= and < sin 2 r/7 osc >= |, and the 
propagation matrix, which describes u a —> Up oscillation, is given by 

Pac/3 — X! I^a«| 2 |^/3i| 2 (3) 

The propagation matrix P a p given by Eq. ([3]) allows to calculate flavour com- 
position at observation for given flavour composition at generation. In particu- 
lar for normal generation composition valid for unsuppressed channels of ir decays 
u e : : v T = 1 : 2 : one obtains from Eq. ([3]) the composition at observation 1:1:1. 



In the work^^ are given many examples of absence equipartition at observations: 
the suppression of muon decays in the chain of pion decay, production of neutrinos 
in neutron beam, decay of neutrinos and others. 



3. Observational effects caused by oscillations 

In this paper we will consider three observational effects caused by oscillations: 
(i) appearance of v e neutrinos important for resonant reaction v e + e~ — > W~ — > 
hadrons, (ii) appearance of r neutrinos and oscillation of mirror to ordinary neutrinos. 
In this section we consider the first two effects; oscillations of mirror neutrinos will 
be discussed in section . 

Resonant interaction of v e neutrinos. 

Glashow considered in 1960 the resonance reaction u e + e~ — > W~ — > pT + zy 
The resonance reaction u e + e~ — > W~ — > hadrons was first suggested in 1977 in^^. 
We will follow here the analytic approach of this paper. 

The resonant production of W~ in v e e collisions occurs at v e energy 

£„ = ^ = 6.3x 10 6 GeV. (4) 
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Integrating over the Breit-Wigner resonance, one obtains analytically^^ 1 the rate of 
resonant events in underground detector with number of electrons N e , given by 

v ves = 2ncr cS E Jp e (E )N e , (5) 

where Jp e (Eo) is the diffuse flux of UHE T> e , 2tt is the solid angle for which deep 
underground detector is open for UHE neutrino flux and a e s is the value left after in- 
tegration over the Breit-Wigner formula, which has meaning of effective cross-section 

Vcs = ^pf F = 2.7 x 10- 32 cm 2 , (6) 

where Gp is the Fermi constant. 

For cosmogenic neutrinos produced in jry collisions with CMB photons, the pro- 
duction of v e is strongly suppressed in comparison with v e neutrinos, because the latter 
are produced in the chain of decay of the positively charged pions: n + — > n + — > u e , 
while v e in the chain of decay of negatively charged pions 7r~ — > fi~ — > v e . Produc- 
tion of 7r + occurs in the resonant reaction p + 7 — >• A + , while production of ir~ goes 
through p + 7— > p + n~ + n~ with a small cross-section. As has been first noted in 
— > v e oscillation considerably increases the flux of cosmogenic v e neutrinos. If 
to take the generation flavour ratio v e : : u T =0:1:0, we obtain from Eq. ([3]) the 
ratio at observation v e : : u T = | : 1 : 1. For this flavour ratio the resonant signal 
is seen well over the background. 



Tau neutrinos 

Tau neutrinos are not produced in the chain of pion decays and appear at observation 
due to oscillations. 

In underground detectors tau-neutrinos produce the characteristic two-bang effect 
the first hadronic shower appears in v T + iV — > r + hadrons interaction, and the 
second one is produced by tau-lepton decay r — > v T + hadrons. At energy of tau- 
neutrino E ~ 10 15 eV the distance between two showers is about 50 m, and such 
event can be observed in IceCube detector. The tau lepton propagating between two 
shower vertex radiates, like muon, the Cherenkov photons and this radiation can be 
also detected. Fig. [T] illustrates the detection of tau- neutrino in the deep-underground 
detector. 
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Figure 1: Double-bang effect produced by tau- neutrino in underground detector 

The UHE neutrinos produce Earth-skimming effect , which can be observed 
by gigantic EAS arrays like the Auger detector. This is a mechanism of detection of 
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Figure 2: Earth-skimming effect produced by UHE tau-neutrinos, which can be observed by gigantic 
EAS detectors 

horizontally moving neutrinos, but tau-neutrinos have the advantages in comparison 
with other neutrinos. Tau-neutrino can cross the large thickness of matter in the 
earth, due to effect of regeneration. Interacting with matter it produces tau-lepton, 
which at large energies propagates to large distance and decays producing again tau- 
neutrino. At energy E > 10 18 eV tau-lepton has a decay length about 50 km and 
it can decay in the air producing hadrons which induce EAS propagating almost 
parallel to the Earth surface (see Fig. [2]). The surface detectors register the EAS 
electromagnetic component from a lower part of EAS. In this way the upper limit on 

'on, 

UHE tau-neutrinos was recently obtained at the Auger detector 1 ^. 



4. Cosmogenic neutrinos in the dip model 



Starting from pioneering work ^ the fluxes of cosmogenic neutrinos have been 
calculated in many works - . The predicted fluxes differ very considerably, de- 
pending on the different assumptions about mass composition of accelerated particles, 
on maximum energy of acceleration and on cosmological evolution of the sources. We 

shall pre sent h ere the UHE neutrino fluxes calculated in the dip model for observed 
UHECR [28l2i 

The dip model for UHECR 

The pair production dip is a feature of interaction of extragalactic UHE protons 
propagating trough CMB radiation. It is caused by energy losses of protons due to 
P+7cmb - ► e + +e~+p scattering. This feature in proton spectrum, in principle, is very 
similar to the GZK cutoff which is caused by photopion production p+7cmb — » N+tt. 
Both features are convenient to analyse in terms of modification factor rj(E) , which 
is defined as a ratio of the diffuse proton spectrum J P (E) calculated with all energy 
losses included to the so-called unmodified spectrum Jp nm (E), when only adiabatic 
energy losses due to expansion of the universe are taken into account: 

The modification factor should be considered as the theoretical spectrum. Being 
defined as a ratio of two spectra it is free from many uncertainties. The modification 
factor is calculated using the power-law generation function Q(E g ) oc E^ 13 , where E g 
is the energy of a proton at generation in a source and 7 ff is the generation index, 
The calculated modification factor gives the excellent agreement with observational 
data for 7 9 = 2.6 — 2.7 and for absence of source evolution, i.e. for evolutionary factor 
(1 + z) m with m = 0. However, inclusion of evolution gives also good agreement with 
the data. This is not a surprising thing, since inclusion of evolution means introducing 
two additional free parameters, m and z max - In Fig. [3] we show the comparison of the 
predicted dip and GZK cutoff with observational data of Akeno-AGASA, Yakutsk, 
HiRes and Auger detectors. 

The observable part of the dip is extended from 1 x 10 18 eV, where modification 
factor reaches 1, and up to 4x 10 19 eV, where the GZK cutoff begins. The agreement of 
the dip with all data is very good. However, one can see that the modification factor 
in the Akeno-AGASA data and in HiRes data exceeds 1 at E < 1 x 10 18 eV. By 
definition the modification factor cannot be larger than 1. The excess of modification 
factor over 1 signals the appearance of the new component of cosmic rays, which can 
be nothing but galactic cosmic rays. Thus , this excess evidences for transition from 
galactic to extragalactic cosmic rays. Strictly speaking at the energy E ~ 1 x 10 18 eV 
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Figure 3: Theoretical pair-production dip and GZK cutoff in comparison with the observational 
data for non-evolutionary models with generation index j g = 2.6 — 2.7. The data of HiRes and 
Auger detectors show steepening of the spectrum consistent with the GZK cutoff. The excess of 
experimental modification factor over ?7 = lat-E<lX 10 18 eV evidences for the new component, 
which is given by galactic cosmic rays. 



transition is completed and at higher energy the cosmic rays are strongly dominated 
by extragalactic component. The transition occurs at lower energy, at the second 
knee seen at the range (4 — 8) x 10 18 eV in the different experiments. The steepening 
observed in the HiRes spectrum is confirmed as the GZK feature by the measured 
value of E1/2 in the integral spectrum. E\ / 2 is a characteristic of the GZK cutoff in the 
integral spectrum ^^21 it j s defined as energy where the integral spectrum calculated 
with all energy losses included becomes half of the power-law extrapolation KE^ 1 
from the low energies. In SSJ an( j CSJ ^ wag demonstrated that E\i2 is a model- 
independent value which equals to 10 19 ' 72 eV. Fig. H] shows how E\/2 was found from 
integral HiRes spectrum. The ratio k = J(> E)/KE^ 1 is plotted as function of 
energy, where J(> E) is measured integral spectrum and KE~"< is its power-law 
extrapolation. The equality k si 1 shows that the chosen approximation is good. The 
energy where integral spectrum crosses the line k = 1/2 gives Ei/ 2 . The corresponding 
value is found Ey2 — 10 19 ' 73±0 07 in excellent agreement with the theoretical value. 
The confirmation of the pair-production dip and the GZK cutoff in observational 
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Figure 4: Ei/ 2 as numerical characteristic of the GZK cutoff in the integral HiRes spectrum. 



data evidences that primary spectrum at E > 1 x 10 18 eV is strongly dominated by 
protons. 

Neutrino fluxes in the dip model 

To calculate the neutrino flux produced by UHE protons it is enough to know 
the generation rate of UHE protons at each cosmological epoch, which we take in 
the form Q(E)(1 + z) m , where factor (1 + z) m describes the cosmological evolution 
of the sources. One should also introduce z max up to which the assumed evolution 
of the sources holds. Q(E) gives the rate of proton generation at z = 0, i.e. the 
number of protons with energy E generated per unit of comoving volume per unit 
time. Expressed in terms of the Lorentz factor T = E/m p and the emissivity £ at 
z = 0, i.e. energy generated per unit comoving volume and unit time, the generation 
rate is given by 

Co 
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where j g = 2.6 — 2.7 is the genera tion in dex and r min ~ 1 is assumed. 

However, following the works l ^ l ^ we assume that acceleration index in each 
source is the same 7 acc = 2.0 — 2.2, but the maximum energies of acceleration E%£? are 
different and distribution of sources over E™ x results in steepening of the generation 
spectrum at E > E c to -y g ~ 2.6 — 2.7. 

The model with the minimum UHE neutrino flux corresponds to absence of the 
evolution m = and low maximum energy of acceleration E™ x = 1 x 10 21 eV. Less 
important is the assumption for value of E c , for which we use E c ~ 1 x 10 18 eV, The 
calculated minimum flux is shown in the upper panel of Fig. [51 This flux is very small, 
but it could be mare inally detectable by IceCube at E > 10 17 eV and by JEM-EUSO 
at E > 10 19 eV. 

One can maximize the flux introducing the cosmological evolution of the sources 



and assuming high maximum acceleration energy. This flux is shown in the lower 
panel of Fig. [5]for the following choice of parameters: m = 4.0, z max = 6.0, 7 9 = 2.45 
and emissivity Cq = 1.2 x 10 46 erg Mpc -1 yr -1 . One may notice the worse agreement 
with the dip with than in the case of the non-evolutionary model. 
Neutrino fluxes from AGN in the dip model 
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Figure 5: UHE neutrino fluxes in the non-evolutionary dip model (upper panel) and in the evolu- 
tionary dip model (lower panel). The neutrino fluxes are shown by curves for the sum of all 
neutrino flavours. The following parameters are used in calculations: m — 0, z max = 2, 7 9 = 
2.7, E max = 1 x 10 21 eV, E c = 1 x 10 18 eV and C = 3.5 x 10 46 erg/Mpc 3 yr (upper panel). In 
the lower panel the neutrino flux is maximized by the following choice of parameters: m = 4.0, 
z m ax = 6.0, 7 9 = 2.45 and emissivity £q = 1-2 x 10 46 erg Mpc~ 1 yr _1 . 



AGN are the most promising sources of the observed UHECR as far as acceleration 
and total energy output is concerned l ^4 | 29 t There are also some obse rvational indi- 
cations in the form of the correlations of UHECR particles with AGN^SESETD yy e 

IOQI 

present here the results of calculations 021 performed in phenomenological approach. 



We assume the generation rate of UHE protons as described above with spectral in- 
dex 7 ff = 2.0 at E < E c and 7 9 = 2.52 above this energy due to assumed distribution 
of AGN over E mauX . The value 7 g = 2.52 is chosen to provide the best fit to the dip 
(see Fig. [6]). The UHE neutrino flux is calculated due to interaction with CMB. The 
evolution of AGN is taken according to X-ray observations (l + z) m with m = 2.7 
up to z c = 1.2. At larger z the evolution is frozen up to z max = 2. Notice, that this 
is very weak evolution, but even such evolution combined with E™ x = 1 x 10 22 eV 
makes UHE neutrino fluxes detectable by EUSO and radio detectors. Neutrino fluxes 
in Fig. [6] are given for one neutrino flavour. 
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Figure 6: UHE neutrino flux in the dip model with AGN as the sources of UHECR. The cosmological 
evolution of AGN with m = 2.7 up to z c = 1.2 is taken from X-ray observations of AGN. At larger 
z the evolution is frozen up to z max = 2.0. The fit of the dip is very good, though requires 7 S = 2.52 
different from the non-evolutionary case m = 0. The neutrino fluxes are given for one neutrino 
flavour and they are detectable especially for the case E"™ x = 1 x 10 22 eV. 



5. Cascade upper limit on UHE neutrino flux 

The Waxman-Bahcall ^3 upper bound is not applicable for UHE neutrinos. This 
bound is obtained from equality of the accompanying UHE proton flux and the ob- 
served flux of UHECR, assuming some relation between proton and neutrino fluxes at 
production. However, one can see from Fig. [5] that UHE neutrino fluxes differ for the 
same proton flux by two-three order of magnitudes at different energies (see upper 
and lower panels). The most interesting case of UHE neutrino flux produced by top- 
down models are not limited by the Waxman-Bahcall bound because accompanying 
proton flux is negligibly small and most produced particles are pions and kaons. 

However, the Waxman-Bahcall limit is useful as a low-flux benchmark for the 
future experiments. 

The most efficient upper bound for UHE neutrinos, applicable for both cosmogenic 



neutrinos and neutrinos from top-down models, is given by the cascade upper limit 
first considered in® (see alsoSD). 

The cascade upper limit on UHE neutrino fluxes I 5 | 42 | 43t j g p rov ided due to e-m 
cascades initiated by UHE photons or electrons which always accompany production 
of UHE neutrinos. Colliding with low-energy target photons, a primary photon or 
electron produce e-m cascade due to reactions 7+7t ar e + + e~, e+7 ta r — > e'+l', etc. 
(see Fig. [7j). The standard case, valid for cosmogenic neutrinos, is given by production 
of UHE neutrinos in extragalactic space, and the cascade develops due to collisions 
with CMB photons (7 tar = 7cmb)- In case neutrino production occurs in a galaxy, the 
accompanying photon can either freely escapes from a galaxy and produce cascade in 
extragalactic space, or produce cascade on the background radiation (e.g. on CMB 
or infrared) inside the galaxy. In the latter case the galaxy should be transparent for 
the cascade photons in the range 10 MeV - 100 GeV. The spectrum of the cascade 
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Figure 7: Electromagnetic cascade developing on the target photons (t), e.g. CMB. 

photons is calculated I 5 | 42 | 43 | 44t . m j QW ener gy p ar ^ ^ j s ^ E~ 3 / 2 , at high energies 
oc E~ 2 with a cutoff at some energy e 7 . The energy of transition between two regimes 
is given approximately by e c ~ (e i /3)(e 7 /m e ) 2 , where e t is the mean energy of the 
target photon. In case the cascade develops in extragalactic space e t = 6.35 x 10 -4 eV, 
e 7 ~ 100 GeV (absorption on optical radiation), and e c ~ 8 MeV. The cascade 
spectrum is very close to the EGRET observations in the range 3 MeV - 100 GeV 
^J. The observed energy density in this range is Regret ~ (2 — 3) x 10 -6 eV/cm 3 . 
It provides the upper limit for the cascade energy density. The upper limit on UHE 
neutrino flux J v (> E) (sum of all flavors) is given by chain of the following inequalities 

47r r 00 47r r°° 4ir 

u cas > — / EJ u {E)dE > —E / J v {E)dE = —EJ U (> E), (9) 

C JE C JE C 

which in terms of the differential neutrino spectrum J U (E) gives 

E 2 J U (E) < —u cas , with u cas < wegret (10) 

47T 

Eq. fllOp gives the rigorous upper limit on the neutrino flux. It is valid for neu- 
trinos produced by HE protons, by topological defects, by annihilation and decays 
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Figure 8: The experimental upper limits on UHE neutrino fluxes in comparison with the e-m cascade 
upper limit in assumption of E~ 2 generation spe ctru m (curve E~ 2 cascade) and with p red ictions 
for cosmogenic neutrinos and neutrinos from TDs — ' . The plot is the modified one from SSI _ 



of superheavy particles, i.e. in all cases when neutrinos are produced through decay 
of pions and kaons. It is valid for production of neutrinos in extragalactic space and 
in galaxies, if they are transparent for the cascade photons. It holds for arbitrary 
neutrino spectrum falling down with energy. If one assumes some specific shape of 
neutrino spectrum, the cascade limit becomes stronger. For E~ 2 generation spectrum, 
which is usually assumed for analysis of observational data one obtains the stronger 
upper limit. Given for one neutrino flavour it reads 

E 2 JAE) < - ~ , ^ -, (11) 

A > ~ 3 47rln(£ max /£ min )' 1 ; 

where i = + or i = v e + v e . 

This upper limit is shown in Fig. [HJ One can see that the observations almost reach 
the cascade upper limit and thus almost enter the region of allowed fluxes. 

The most interesting energy range in Fig. [8] corresponds to E v > 10 21 eV, where 
acceleration cannot provide protons with sufficient energy for production of these neu- 
trinos. In fact this statement is valid only for shock acceleration. In principle, e.g. in 
AGN the different mechanisms of acceleration might operate, such as unipolar induc- 
tion or pinch-instability, and they can provide the higher energy of acceleration, but 
these mechanisms are not developed enough for numerical calculations of produced 
fluxes of accelerated particles. At present the region of E v > 10 21 eV, and especially 
E v ^> 10 21 eV is considered as a signature of top-down models, which provide these 
energies quite naturally. Below we shall consider three top-down models: Superheavy 



Dark Matter, Topological Defects and Mirror Matter. 



6. UHE neutrinos from Superheavy Dark Matter (SHDM) 

SHDM is one of the models for cosmological cold dark matter I 47 | 48 t The mQ st 
attractive mechanism of production is given by creati on of superheavy particles in 
time- varying gravitational field in post-inflation epoch I 49 | 50t Creation occurs when 
the Hubble parameter is of order of particle mass H(t) ~ mj. Since the maximum 
value of the Hubble parameter is limited by the mass of the inflaton H(t) < ~ 
10 13 GeV, the mass of X-particle is limited by m^, too. For example, mx ~ 3 x 
10 13 GeV results in Qxh 2 ~ 0.1, as required by WMAP measurements. 

Being protected by some symmetry, SHDM particles with such masses can be 
stable or quasi-stable. In case of gauge symmetry they are stable, in case of gauge 
discrete symmetry they can be stable or quasi-stable. Decay can be provided by 
superweak effects: wormholes, instantons, high- dimension operators etc. 

Like any other form of cold dark matter, X-particles are accumulated in the halo 
with overdensity 2.1 x 10 5 . 

SHDM particles can produce UHECR and high energy neutrinos at the decay 
of X-particles (when the protecting symmetry is broken) and at their annihilation, 
w hen the sy mmetry is exact. The scenario with decaying X-particles was first studied 
in I 47 | 51 | 52 t ^ n interesting scenario with stable X-particles, when UHE particles 
are produced by annihilation of X-particles has been put forward in ^3}. In this 
scenario superheavy X-particles have the gauge charge and they are produced at post- 
inflationary epoch by close pairs, forming the bound systems. Loosing the angular 
momentum, these particles inevitably annihilate in a close pair. 

The UHE particles (protons, pions and neutrinos from the chain of pion decays) 
are produced as a result of QCD cascading of partons. The calculations of fluxes and 
spectra are nowadays reliably performed by Monte Carlo ^* and using the DGLAP 
equations - i^u. The spectra of protons, photons and neutrinos are shown in 
Fig. [9] for the case of SHDM particles with mass Mx = 1 x 10 14 GeV. The spectrum 
of photons is normalized by the AG AS A excess. In case it is absent, like in HiRes 
and Auger observations, all fluxes, including neutrino, must be lowered by factor 3 - 
5. 



7. UHE neutrinos from Topological Defects (TDs) 

As has been first noticed by D. A. Kirzhnitz each spontaneous symmetry 
breaking in the early universe is accompanied by the phase transition. Like the phase 
transitions in liquids and solids, the cosmological phase transitions can give rise to 
topological defects (TDs), which can be in the form of surfaces (cosmic textures), lines 
(cosmic strings) and points (monopoles). In many cases TDs become unstable and 
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Figure 9: Spectra of neutrinos (upper curve), photons (middle curve) and protons (two lower curves) 
in SHDM model compared with AGASA data, according to calculations of . The neutrino 
flux is dominated by the halo component with small admixture of extragalactic flux. The flux of 
extragalactic protons is shown by the lower curve (extr. p). The fluxes are normalized by the 
AGASA excess at E > 1 x 10 20 eV. In case it is absent, like in HiRes and Auger data, all fluxes, 
including neutrino, must be lowered by factor 3-5. 



decompose to constituent fields, superheavy gauge and Higgs bosons (X-particles), 
which then decay producing UHECR. It could happen, for example, when two seg- 
ments of ordinary string, or monopole and antimonopole touch each other, when 
electrical current in superconducting string reaches the critical value and in some 
other cases. The decays of these particles, if they heavy enough, produce particles of 
ultrahigh energies including neutrinos. 

The following TDs are of interest for UHECR and neutrinos : 
monopoles (G — > H X U(l) symmetry breaking), ordinary strings {U(l) symmetry 
breaking) with important subclass of superconducting strings, monopoles connected 
by strings (G —>■ H x U(l) symmetry breaking with subsequent U(l) — > Zn symmetry 
breaking, where Zn is discrete symmetry). The important subclass of the monopole- 
string network is given by necklaces, when Z^ = Z 2 , i.e. each monopole is attached 
to two strings. We shall shortly describe the production of UHE particles by these 
TDs. (see^" for more details). 

(i) Ordinary strings. 

There are several mechanisms by which ordinary strings can produce UHE particles. 
For a special choice of initial conditions, an ordinary string loop can collapse to 
a double line, releasing its total energy in the form of X-particles. However, the 
probability of this mode of collapse is extremely small, and its contribution to the 



overall flux of UHE particles is negligible. 

String loops can also produce X-particles when they self- intersect. Each intersec- 
tion, however, gives only a few particles, and the corresponding flux is very small. 

The loops undergo oscillation and an important property of it is the periodic 
appearance of the cusps, the loop points with velocity of light. The near-cusp seg- 
ments, moving with large Lorentz factor may overlap and annihilate producing the 
constituent superheavy particles with large Lorentz factors. Therefore the particles 
from X-decays are further boosted by large Lorentz factor. The energy released in a 
single cusp event can be quite large, but again, the resulting flux of UHE particles is 
too small to account for the UHECR observations, 
(ii) Superconducting strings. 



As was first noted by WitterF*, in a wide class of elementary particle models, strings 
behave like superconducting wires. Moving through cosmic magnetic fields, such 
strings develop electric currents. Superconducting strings produce X-particles when 
the electric current in the strings reaches the critical value. This process is strongly 
increased in cusps where in a small fraction of loop the current becomes supercritical 
and X-particles, the charge carriers, leave the string and are decaying to the high- 
energy ordinary particles. Their energies are further boosted by cusp Lorentz factor. 
Superconducting strings cannot be the sources of observed UHECR ^SJ because of 
the absorption of the particles on the way from a source to observer, but they can 
produce the observable flux of UHE neutrinos. 

(in) Network of monopoles connected by strings. 
The sequence of phase transitions 



results in the formation of monopole-string networks in which each monopole is at- 
tached to N strings. Most of the monopoles and most of the strings belong to one 
infinite network. The evolution of networks is expected to be scale-invariant with 
a characteristic distance between monopoles d = Kt, where t is th e ag e of Universe 
and k = const. The production of UHE particles are considered in 1 ^. Each string 
attached to a monopole pulls it with a force equal to the string tension, \i ~ 
where n s is the symmetry breaking vev of strings. Then monopoles have a typical 
acceleration a ~ /x/m, energy E ~ [id and Lorentz factor T m ~ fid/m, where m is the 
mass of the monopole. Monopole moving with acceleration can, in principle, radiate 
gauge quanta, such as photons, gluons and weak gauge bosons, if the mass of gauge 
quantum (or the virtuality Q 2 in the case of gluon) is smaller than the monopole 
acceleration. The typical energy of radiated quanta in this case is e ~ T m a. This 
energy can be much higher than what is observed in UHECR. However, the produced 
flux (see^J 1 ) is much smaller than the observed one. 
(vi) Necklaces. 

Necklaces are hybrid TDs corresponding to the case N = 2 , i.e. to the case when 




H x 17(1) ^ H x Z 1 



(12) 



each monopole is attached to two strings. This system resembles "ordinary" cosmic 
strings, except the strings look like necklaces with monopoles playing the role of 
beads. The evolution of necklaces depends strongly on the parameter 



r = m/fj,d, 



(13) 



where m is a mass of a monopole, [i is mass per unit length of a string (tension of a 
string) and d is the average separation between monopoles and antimonopoles along 
the strings. As it is argued in Ref. necklaces might evolve to configurations with 
r ^> 1. Monopoles and antimonopoles trapped in the necklaces inevitably annihilate 
in the end, producing first the heavy Higgs and gauge bosons (X-particles) and then 
hadrons. The rate of X-particle production can be estimated as ^~ 



2 

r • ji 
t 3 mx 



(14) 



This rate determines the rates of pion and neutrino production with energy spectrum 
calculated in Ref. . Restriction due to e-m cascade radiation demands the cascade 
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Figure 10: Diffuse spectra of neutrinos, protons and photons from necklaces. The upper curve 
shows neutrino flux, the middle - proton flux and two lower curves - photon fluxes for two cases of 
absorption. The thick curve gives the sum of the proton and the higher photon flux. 



energy density u cas < 2 • 10 6 eV/cm 3 . The cascade energy density produced by 
necklaces can be calculated as 



UJ, 



1 f 2 f 
2 Ur No 



*o dt 



1 



t 3 (1 + z) 



3 

I" 



aJk 1 .21 



(15) 



where f w ~ 0.5 is a fraction of total energy release transferred to the cascade. There- 
fore, r 2 fi and the rate of X-particle production (TH|) is limited by cascade radiation. 

The fluxes of UHE protons, photons and neutrinos from necklaces are shown in 
Fig. [10] according to calculations of The mass of X-particle is taken mx = 
1 x 10 14 GeV. Neutrino flux is noticeably higher than in the case of conservative 
scenarios for cosmogenic neutrinos and neutrinos from SHDM. 

In the recent work^l n was indicated that in some models of necklaces there could 
be fast annihilation of monopoles and separation of pairs monopole-antimonopole con- 
nected by strings from the necklaces. In these models UHE neutrino flux is suppressed. 
However, in some other models the detectable neutrino flux can exist. 



8. Mirror matter and mirror neutrinos 



Existence of mirror matter in our universe is the old idea which was put forward 
in the end of 1950s. Mirror matter can be most powerful source of UHE neutrinos not 
limited by e-m cascade limit. Produced as the mirror neutrinos, they can oscillate 
into ordinary neutrinos. All mirror particles that accompany production of mirror 
neutrinos in the mirror matter remain invisible for our detectors. 

Concept of mirror matter 

Mirror matter is based on the theoretical concept of the space reflection, as first 
suggested by Lee and Yangl^ in 1956 and developed by Landau in 1957, Salam 

in 1957, and most notably by Kobzarev, Okun, Pomeranchuk in 1966 (see 
recent exciting historical review by Okun^^). 

This concept can be explained in the following way. 

The Hilbert particle space is assumed to be a representation of the extended 
Lorentz group, which includes the space coordinate reflection x — > —x. Since the 
coordinate operations, reflection x — > — x and time shift t — ► t + At, commute, the 
corresponding operations in the particle space, I r and Hamiltonian H, must commute, 
too: 

[H,I r ] = (16) 

It implies that operator I r must correspond to the conserved value. Since parity 
P according to assumption of Lee and Yang is not conserved, I r should be defined 
somehow else. Lee and Yang suggested I r = P x R, where R transfers particle to 
mirror particle. Since parity operator P interchange left and right states, one obtains 

I r *L = K and I r ^ R = ^' L , (17) 

where primes indicate the states in mirror particle space. 

The assumption of Landau^ 1 was R = C, i.e. one may say that he suggested to 
use antiparticles as the mirror space, but then CP must be conserved which as we 



know today is not the case. 

Oscillation of mirror and ordinary neutrinos 

Kobzarev, Okun, and Pomeranchuk, ^SJ suggested that ordinary and mirror sectors 
communicate only gravitationally. For the description of this interaction one can use 
dimension 5 operator obtained as SU(2)l X !7(1) x SU{2)' R x U{1)' scalar: 

£comm = T-^K^X^'); (18) 

where Mp\ = 1.2 x 10 19 GeV is the Planckian mass, which implies the gravitational 
interaction, and <ft, <$' are the electroweak Higgses from visible and mirror sectors, 

After spontaneous electroweak symmetry breaking (0) = (</>') = v the Lagrangian 
( {TBI) generates the terms, which mix visible and sterile neutrinos, 

v 2 

£ mix = TT~ vv , ( 19 ) 
Mpi 

where v— 174 GeV is vacuum expectation value of Higgses, and /x = v 2 /Mp\ = 
2.5 x 10 -6 eV is the mixing parameter. 

Eq. (1191) implies oscillation of mirror and ordinary neutrinos. This oscillation 
described by Eq. ffT9]) . but with M not necessarily being the Planckian mass, was 
first suggested by Berezhiani and Mohapatra and Foot and Volkas^^ 1 in 1995. 

UHE neutrinos from mirror TDs 

Any cosmological scenario for mixed ordinary and mirror matter must provide the 
suppression of the mirror matter and in particular the density of mirror photons and 
neutrinos at the epoch of nucleosynthesis. It can be obtained in the two-inflaton model 
The rolling of two inlatons to minimum of the potential is not synchronized, and 
when the mirror inflaton reaches minimum, the ordinary inflaton continues its rolling, 
inflating thus the mirror matter produced by the mirror inflaton. While mirror matter 
density is suppressed, the mirror topological defects in two-inflatons scenario with 
curvature- driven phase transition can strongly dominate Mirror TDs copiously 
produce mirror neutrinos with extremely high energies typical for TDs, and they are 
not accompanied by any visible particles. Therefore, the upper limits on HE mirror 
neutrinos in our world do not exist. All HE mirror particles produced by mirror 
TDs are sterile for us, interacting with ordinary matter only gravitationally, and only 
mirror neutrinos can be efficiently converted into ordinary ones due to oscillations. 
The oscillations of mirror and visible neutrinos in the gravitational-mixing scenario 
has been studied in detail in^D. The probability of oscillation of mirror neutrino u' 
into visible neutrinos is large. In particular, for v' neutrino it is given by 

P^ Ve = ~ sin 2 29 12 , P V ^=P V>T = 1 - 1 sin 2 2# 12 , Y. P ^ = \ ( 20 ) 



9. Conclusions 



UHE neutrino astronomy is characterised by well balanced program. 

As the secure part of this program, there are cosmogenic neutrinos, for production 
of which the particle beam (UHECR particles) and the target (CMB photons) are 
well known. In the case when UHECR primaries are protons and the observed dip 
is a feature of UHE proton interaction with CMB, there is the robustly predicted 
lower limit for UHE neutrino flux (see upper panel of Fig. EJ), which is marginally 
detectable by IceCube at E < 1 x 10 18 eV and by JEM-EUSO at E > 1 x 10 19 eV. This 
lower limit is supported by the following observations: the features of UHE proton 
interaction with CMB, GZK cutoff and dip, are found, mostly in HiRes observations, 
and the fluorescent data of HiRes evidence for proton-dominated mass composition. 
However, the Auger fluorescent data favour the mixed nuclei composition. In this 
case the cosmogenic neutrinos can be undetectable for existing projects. 

The large fluxes of cosmogenic neutrinos, detectable by existing projects, cor- 
respond (at the fixed UHECR flux) to cosmological evolution of the sources, flat 
generation spectra and large maximum energy of acceleration. The lower panel of 
Fig. [5] presents UHE neutrino fluxes for the extreme hypothetical assumptions: very 
large E max and strong evolution of the sources up to z max = 6. 

The fundamental problem of astrophysics involved in prediction of UHE cos- 
mogenic neutrinos is acceleration of particles. The shock acceleration at present 
knowledge of its theory cannot provide maximum energy of acceleration higher than 
10 21 — 10 22 eV, and thus energies of cosmogenic neutrinos do not exceed 3 x 10 20 eV. 
However, in case of AGN as UHECR sources, the alternative mechanisms of acceler- 
ation, such as unipolar induction and plasma mechanisms, can work and accelerate 
particles to much higher energies. Unfortunately, these mechanisms are not devel- 
oped enough for numerical calculations of produced rates and spectra at generation. 
The detection of UHE neutrinos from the individual sources can help to solve this 
fundamental problem of astrophysics. 

The second part of UHE neutrino astronomy is provided by top-down sources: 
SHDM and topological defects in ordinary and mirror matter. The decays of unstable 
constituent fields of TDs can produce neutrinos with extremely high energies up to 
GUT scale and above. This is the general property of top-down models. However, 
the fluxes are very model dependent, and in case of TDs they differ much for different 
types of topological defects. There could be no flux estimates based on the general 
common properties of all TDs: some of them give low undetectable fluxes, from some 
exceptional TDs one can expect the detectable fluxes. For long time the necklaces 
have been considered as most promising TDs for UHECR and neutrinos^D 1 , but recent 
work ^3 p U t s the doubts on these topological defects as sources of UHE particles in 
recent cosmological epochs. The signature of top-down models is very high neutrino 



energies, unreachable for accelerator neutrinos. These neutrinos are most reliably 
detectable by radio and acoustic methods, as well as by EUSO detectors. 

In both cases of cosmogenic and top-down models, the neutrino fluxes are con- 
strained by the cascade upper limit. As exception, the mirror neutrinos do not respect 
this limit, and their fluxes can be even larger. 

The search for UHE neutrinos in any case is a search for a new physics, either 
for astrophysics (the new acceleration mechanisms and cosmological evolution of the 
sources) or for topological defects, mirror topological defects and superheavy dark 
matter. 
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